Determination of the pion charge form factor for (3 2 =0.60-1.60 GeV 2 
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The data analysis for the reaction 1 H(e, e'7r + )n, which was used to determine values for the charged 
pion form factor for values of Q 2 =0.6-1.6 GeV 2 , has been repeated with careful inspection of all 
steps and special attention to systematic uncertainties. Also the method used to extract FV from 
the measured longitudinal cross section was critically reconsidered. Final values for the separated 
longitudinal and transverse cross sections and the extracted values of F„ are presented. 

PACS numbers: 14.40.Aq,11.55.Jy,13.40.Gp,25.30.Rw 



I. INTRODUCTION 



Hadron form factors are an important source of information on hadronic structure. Of these, the electric form 
factor, F n , of the charged pion plays a special role. One of the reasons is that the valence structure of the pion is 
relatively simple. The hard part of the ir + form factor can be calculated within the framework of perturbative QCD 
(pQCD) as the sum of logarithms of Q 2 multiplied by powers of l/Q 2 As Q 2 — > oo, only the leading-order term 
remains 

where a s is the strong-coupling constant and f n the pion decay constant. Thus, in contrast to the nucleon case, the 
asymptotic normalization of the pion function is known from the decay of the pion. The theoretical prediction for 
^ at experimentally accessible Q 2 is less certain, since the calculation of the soft contributions is difficult and model 
dependent. This is where considerable theoretical effort has been expended in recent years. Some examples include 
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next-to-leading order (NLO) QCD @,J1, QCD Sum Rules 0,(1], Constituent Quark Models @, and Bethe-Salpeter 
Equation [7J calculations. (See Ref. [8( for a review.) Some of these approaches are more model independent than 
others, but it is fair to say that all benefit from comparison to high quality F^ data, to delineate the role of hard 
versus soft contributions at intermediate Q 2 . 

The experimental measurement of the pion form factor is quite challenging. At low Q 2 , can be measured in a 
model independent manner via elastic scattering of n + from atomic electrons such as at the CERN SPS ||. Above 
Q 2 > 0.3 GeV 2 , one must determine F v from pion electroproduction on the proton. The dependence on F T enters the 
cross section via the t-channel process, in which the incident electron scatters from a virtual pion, bringing it on-shell. 
This process dominates near the pion pole at t = m 2 . The physical region for t in pion electroproduction is negative, 
so measurements should be performed at the smallest attainable values of — t. To minimize background contributions, 
it is also necessary to separate the longitudinal cross section a^, via a Rosenbluth L/T(/LT/TT) separation. The 
value of F- K {Q 2 ) is then determined by comparing the measured longitudinal cross section at small values of — t, where 
it is dominated by the t-pole term, which contains F v , to the best available model for the reaction 1 H(e,e'ir + )n, 
adjusting the value of F^ in the latter. 

Using the electroproduction technique, the pion form factor was studied for Q 2 values from 0.4 to 9.8 GeV 2 at 
CEA/Cornell [10] and for Q 2 =0.35 and 0.70 GeV 2 at DESY [UG3. Ref. [3 performed a longitudinal/transverse 
separation by taking data at two values of the electron energy. In the experiments done at CEA/Cornell, this was 
done in a few cases only, but even then the resulting uncertainties in <tl were so large that the L/T separated data 
were not used, and ctl was determined by assuming a certain parametrization for a^. Consequently, the values of F^ 
extracted from these data have sizable systematic uncertainties. 

More recently, the 1 H(e, e'ir + )n reaction was measured at the Thomas Jefferson National Accelerator Facility (JLab) 
in order to study the pion form factor from Q 2 =0.6-1.6 GeV 2 . Because of the excellent properties of the electron beam 
and experimental setup at JLab, L/T separated cross sections were determined with high accuracy. These data were 
used to determine the value of F^ and the results were published in Ref. [13j . Since then, the whole analysis chain has 
been repeated with careful investigation of all steps, including the contribution of various systematic uncertainties to 
the final uncertainty of the separated cross sections. Furthermore, the method to determine F„ from the longitudinal 
cross sections was re- investigated, leading to a different method to extract F„. In this paper, we report on these 
studies and present final results for the longitudinal and transverse cross sections, as well as the extracted values of 
F^ from these data. We also discuss in detail the extraction of F n from the measured cross sections, and the related 
extraction uncertainties (model dependence). 



II. EXPERIMENT AND CROSS SECTION DATA ANALYSIS 



The cross section for pion electroproduction can be written as 



rvJ(t)^, (2) 



dE'dtle'd^ ^ ' dtd<p'' 

where 

r 1 w 2 -m 2 

2n 2 E e Q 2 \-e 2M [ ' 

is the virtual photon flux factor, (f> is the azimuthal angle of the outgoing pion with respect to the electron scattering 
plane, t = (p n — q) 2 is the Mandelstam variable, J is the Jacobian for the transformation from dO w to dtd(f> and W 
is the photon- nucleon invariant mass. 

The two-fold differential cross section can be written as 

d 2 a dcr L da T r— dcr LT 

2tt — — ■ = e — h\/2e e+1 — - — cosm 

dtdcf) dt dt v v ' dt r 

+e—^cos2<f>, (4) 
where e = ( I + 2yU- tan 2 % J is the virtual-photon polarization parameter. The ax = cross sections depend 



on W, Q 2 and t. By using the (^-acceptance of the experiment and taking data for the same (central) kinematics 
{W,Q 2 ,t) at two energies, and thus two values of e, the cross sections ctl, <Xt 5 0lt and ott can all be determined. 

Using 2.4-4 GeV electron beams impinging upon a liquid hydrogen target, data for the 1 H(e, e'ir + )n reaction were 
taken at a central value of W — 1.95 GeV for central Q 2 — values of 0.6, 0.75, 1.0 and 1.6 GeV 2 . The scattered electron 
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FIG. 1: (Color online) Representative missing mass distribution. Cuts at 0.925 and 0.960 GeV were applied to select the recoil 
neutron final state. 



was detected in the Short Orbit Spectrometer (SOS) and the produced pion in the High Momentum Spectrometer 
(HMS) of Hall C. 

The data analysis is an updated version of that in Ref. (l3j ]. First, experimental yields were determined. Electrons 
in the SOS were identified by using the combination of a lead glass calorimeter and gas Cerenkov detector. Pion 
identification in the HMS was accomplished by requiring no signal in a gas Cerenkov detector and by using time of 
flight between two scintillator hodoscope planes. The momenta of the scattered electron and the pion at the target 
vertex were reconstructed from the wire chamber information of the spectrometers, correcting for energy loss in the 
target. From these, the values of Q 2 , W, t, and the missing mass were reconstructed. A cut on the latter of 0.925 
to 0.96 GeV was used to select the neutron exclusive final state, excluding additional pion production (Fig. [T]). 
Experimental yields as function of Q 2 , W, t and <fi were then determined by subtracting random coincidences (varying 
with bin but typically 1.2%) and aluminum target window contributions (typically 0.6%) and correcting for trigger, 
tracking and particle- identification efficiencies, pion absorption, local target-density reduction due to beam heating, 
and dead times. Details of these procedures are similar to those found in Ref. [i~4j . 

Cross sections were obtained from the yields using a detailed Monte Carlo (MC) simulation of the experiment, 
which included the magnets, apertures, detector geometries, realistic wire chamber resolutions, multiple scattering 
in all materials, optical matrix elements to reconstruct the particle momenta at the target from the information of 
the wire-chambers of the spectrometers, pion decay (including misidentification of the decay muon as a pion), and 
internal and external radiative processes. 

Calibrations with the over-determined 1 H(e,e'p) reaction were instrumental in various applications. The beam 
momentum and the spectrometer central momenta were determined absolutely to better than 0.1%, while the incident 
beam angle and spectrometer central angles were determined with an absolute accuracy of about 0.5 mrad. The 
spectrometer acceptances were checked by comparison of the data to the MC simulations. Finally, the overall absolute 
cross section normalization was checked. The calculated yields for e + p elastics agreed to better than 2% with 
predictions based on a parameterization of the world data [15[ . 

In the pion production reaction, the experimental acceptances in W , Q 2 and t are correlated. By using a realistic 
cross section model in the MC simulation, possible errors resulting from averaging the measured yields when calculating 
cross sections at average values of W, Q 2 and t, can be minimized. A phenomenological cross section model was 
obtained (see below) by fitting the different cross sections ax of Eqn. (01 globally to the data as a function of Q 2 and 
t in the whole range of Q 2 . The dependence of the cross section on W was assumed to follow the phase space factor 
(W 2 — M 2 )~ 2 , which is supported by previous data 

The experimental cross sections can then be calculated from the measured and simulated yields via the relation 

da{W,Q 2 ,t,cp) \ _ (Y exp ) / da(W,Q 2 ,t,(f>) 
dt I (Y MC ) \ dt 

/ exp \ 

This was done for five bins in t at each of the four Q 2 -values. Here, (Y) indicates that the yields were averaged over 

the W and Q 2 acceptance, W and Q being the acceptance (of high and low e together) weighted average values for 
that i-bin. By using these average values, possible errors due to extrapolating the MC model cross section used to 
outside the region of the experimental data, is avoided. 

By combining for every i-bin (and for the four values of Q 2 ) the </>-dependent cross sections measured at two values 
of the incoming electron energy, and thus of e, the experimental values of ox, 0t, <tlt and cttt can be determined 
by fitting the <j> and e-dependence (Fig. [2|). In this fit, the leading order sin 9 (sin 2 0) of olt(ott), where 9 is the 
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FIG. 2: Example of the measured cross sections, as a function of <fi at Q 2 =0.6 GeV 2 for two values of e. The curves 

shown represent the model cross section used in the Monte Carlo simulation. 
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FIG. 3: (Color online) Separated cross sections ol [solid] and ot [open]. The error bars represent the combination of statistical 
and t uncorrelated systematic uncertainties. In addition, there is an overall systematic uncertainty of about 6%, mainly from 
the t correlated, e uncorrelated systematic uncertainty. The solid and dashed curves denote VGL model calculations for ox 
and or with parameters A 2 =0.393, 0.373, 0.412, 0.458 GeV 2 for Q 2 =0.6-1.6 GeV 2 , and with common A 2 =1.5 GeV 2 . The 
discontinuities in the curves result from the different average W and Q 2 of each i-bin. 



angle between the three- momentum transfer and the direction of the outgoing pion, was taken into account. 1 Those 
values were then used to improve upon the model cross section used in the MC simulation. This whole procedure was 
iterated until the values of ox, ax, <J\jt and cttt converged. The dependence of ax (and a?) on the MC input model 
was small (see below). 

The separated cross sections ax and ax are shown in Fig. [3l They are presented as differential cross sections da/dt 
as a function of t, at the center of the i-bin. The longitudinal cross section exhibits the expected i-pole behavior. The 
transverse cross section is mostly flat. 

The total uncertainty in the experimental cross sections is a combination of statistical and systematic uncertainties. 
All contributions to the systematic uncertainty were carefully investigated, also using the results of extensive single- 
arm L/T separation experiments and of 1 Yi{e,e'p) calibration reactions in Hall C [la ]. The experimental systematic 



1 In the previous analysis [T3| . first oxt and <ttt were determined by adjusting their values (plus a constant term) until the ratios 
Y exp /YMC were constant as function of 6 and <j>. After that, ox and <tt were determined in a Rosenbluth separation. The present 
method is more straightforward and has the advantage that the uncertainties in the separated cross sections are obtained more directly. 
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FIG. 4: Differential cross section comparison between our earlier Q 2 =0.75 GeV 2 results [open circles] (Ref. [13|]) and this work 
[filled circles]. The unseparated cross sections (a uns = eaL + ctt) at high and low e are nearly identical, but the differences 
between the separated ol and ot are somewhat larger. The a un3 error bars include the statistical and epsilon uncorrelated 
systematics only and are in many cases smaller than the plotting symbols. Those for ctl,t include the contribution of all 
statistical and systematic uncertainties. 



uncertainties include contributions that, like the statistical uncertainties, are uncorrelated between the measurements 
at the two e values, and others that are correlated. Most of the uncorrelated ones are common to all t bins, but 
there is a small contribution, estimated as 0.7%, that is also uncorrelated in t. The e-uncorrelated uncertainties in 
<7l are inflated by the factor 1/Ae in the L/T separation, where Ae is the difference (typically 0.3) in the photon 
polarization between the two measurements. The effect on <tt depends on the exact e values. The e-uncorrelated 
systematic uncertainty in the unseparated cross sections common for all i-bins was estimated to be 1.7%, while the 
total correlated uncertainty is 2.8-4.1%, depending on t. Apart from a dependence of the separated cross sections on 
the MC model used, which ranges from 0.2% to a maximum of 3% for one highest t-bin, the largest contributions 
are: the detection volume (1.5%), dependence of the extracted cross sections on the momentum and angle calibration 
(1%), target density (1%), pion absorption (1.5%), pion decay (1%), the simulation of radiative processes (1.5%), and 
detector efficiency corrections (1%). The overall uncertainty is slightly smaller than used in Ref. fl3| . 

The unseparated cross sections and hence also the values of ctl and ctt of the present analysis differ from those of 
our earlier analysis presented in Ref. [l3j . Compared to that analysis, small adjustments were made in the values of 
cuts and efficiencies. Also, a small mistake was found in calculating the value of 9, which affects the calculation of 
the cross section in Eqn. \5\ Finally, as mentioned, the method to separate the cross sections was changed. The cross 
sections in Table fl] are our final values. Except for a few cases, the difference with the previous values is well within 
the total uncertainty quoted in Ref. [13]. As an example, the old and the new cross sections for the case Q 2 = 0.75 
GeV 2 are shown in Fig. [4] It can be seen that the differences in the extracted unseparated cross sections (top panels) 
are very small, but the L/T separation increases them. On average over the four Q 2 cases, <tl is 6% smaller than in 
Ref. [13| and ctt is 3% larger. The largest differences occur for Q 2 =1.0 GeV 2 , where ctl is 14% smaller and or is 
10% larger. 



III. EXTRACTION OF F^(Q 2 ) FROM THE DATA 

It should be clear that the differential cross sections ctl versus t over some range of Q 2 and W are the actual 
observables measured by the experiment. The extraction of the pion form factor from these cross sections can be done 
in a number of approaches, each with their own merits and associated uncertainties. 

Frazer [l7| originally proposed that F n be extracted from <tl via a kinematic extrapolation to the pion pole, and that 
this be done in an analytical manner, a la Chew-Low (l8| . This extrapolation procedure fails to produce a reliable 
answer, since different polynomial fits, each of which are equally likely in the physical region, differ considerably 
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Q 2 W -t 

(GoV 2 ) (GoV) (GoV 2 ) 


dtJL/dt 

((jb/GoV 2 ) 


dor/di 

( fl b/GoV 2 ) 


0.526 1.983 0.026 
0.576 1.956 0.038 
0.612 1.942 0.050 
0.631 1.934 0.062 
0.646 1.929 0.074 


31.360 ± 1.602, 1.927 
24.410 ± 1.119, 1.774 
20.240 ± 1.044, 1.583 
14.870 ± 1.155, 1.366 
11.230 ± 1.469, 1.210 


8.672 ± 1.241 
10.660 ± 1.081 
10.520 ± 1.000 
10.820 ± 0.992 
10.770 ± 1.097 


0.660 1.992 0.037 
0.707 1.961 0.051 
0.753 1.943 0.065 
0.781 1.930 0.079 
0.794 1.926 0.093 


20.600 ± 1.976, 1.895 
16.280 ± 1.509, 1.788 
14.990 ± 1.270, 1.573 
11.170 ± 1.214, 1.416 
9.949 ± 1.376, 1.277 


9.812 ± 1.532 
10.440 ± 1.344 
8.580 ± 1.150 
9.084 ± 1.091 
8.267 ± 1.110 


0.877 1.999 0.060 
0.945 1.970 0.080 
1.010 1.943 0.100 
1.050 1.926 0.120 
1.067 1.921 0.140 


14.280 ± 1.157, 1.103 
11.840 ± 0.887, 0.978 
9.732 ± 0.773, 0.837 
7.116 ± 0.789, 0.747 
4.207 ± 1.012, 0.612 


7.084 ± 0.791 
6.526 ± 0.657 
5.656 ± 0.572 
5.926 ± 0.570 
5.802 ± 0.656 


1.455 2.001 0.135 
1.532 1.975 0.165 
1.610 1.944 0.195 
1.664 1.924 0.225 
1.702 1.911 0.255 


5.618 ± 0.431, 0.442 
4.378 ± 0.356, 0.390 
3.191 ± 0.322, 0.351 
2.357 ± 0.313, 0.310 
2.563 ± 0.356, 0.268 


3.613 ± 0.294 
3.507 ± 0.257 
3.528 ± 0.241 
3.354 ± 0.228 
2.542 ± 0.227 



TABLE I: Separated cross sections ox and <tt from this work. The two listed uncertainties for ox are the combination of 
statistical and f-uncorrelated systematic uncertainties, and the combination of the e-correlated and uncorrelated uncertainties. 
The statistical and t-uncorrelated uncertainties are applied before fitting the VGL model to the data, while the e-correlated 
and uncorrelated uncertainties are applied afterwards. The listed errors for or include all experimental uncertainties. 



when continued to t = m 2 . Some attempts were made to reduce this uncertainty by providing some theoretical 
constraints on the behavior of the pion form factor in the unphysical region, but none proved adequate. 

Bebek et al. [1(3] embraced the use of theoretical input when they used the Born term model of Berends to 
perform a form factor determination. Brauel et al. [l2j similarly used the Born term model of Gutbrod and Kramer 
m| to extract F w . The presence of the nucleon and its structure complicates the theoretical model used, and so an 
unavoidable implication of this method is that the extraction of the pion form factor becomes model dependent. 

As in Ref. [H, H2], the Regge model by Vanderhaeghen, Guidal and Laget (VGL, Ref. [23]) is used here to extract 
Fir- In this model, the pole-like propagators of Born term models are replaced with Regge propagators, and so the 
interaction is effectively described by the exchange of a family of particles with the same quantum numbers instead 
of the exchange of one particle. The model was first applied to pion photoproduction. Most of the model's free 
parameters were determined from data on nucleon resonances. The use of Regge propagators, and the fact that 
both the 7T (J = 0) and the p (J = 1) trajectories are incorporated in the model proved to be essential to obtain 
a good description of the W- and f-dependence of the photoproduction data for both ir + and ir~ particles. For 
electroproduction, the pion form factor and the pirj form factor are added as adjustable parameters, parameterized 
with a monopole form 

F 7V (Q 2 ) = [l + Q 2 /Al}-\ (6) 

The Regge model does a superior job of describing the t dependence of the differential pion electroproduction cross 
sections of [TTI . [l2| than the Born term model. Over the range of —t covered by this work, ctl is completely determined 
by the ir trajectory, whereas o"t is also sensitive to the p exchange contribution. The value of A 2 is poorly known, 
while A 2 is much better known and in the end is determined by the fitting of the model to the o"l data. 

The VGL model for certain choices of A 2 , and A 2 is compared to our data in Figure [3] The VGL cross sections have 

been evaluated at the same W and Q 2 values as the data, resulting in the discontinuities shown. The model strongly 
underestimates or for any value of A 2 used (variation of A 2 within reasonable values can change o"t by up to 40%). 
Since the JLab data have been taken at relatively low values of W « 1.95 GeV, this may be due to contributions 
from resonances, enhancing the strength in or- No such terms are included explicitly in the Regge model. The VGL 
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model calculation for ol gives the right magnitude, but the t dependence of the data is somewhat steeper than that 
of the calculations. This is most visible at Q 2 =0.6 GeV 2 . As in the case of <tt, the discrepancy between the data and 
VGL is attributed to resonance contributions. This is supported by the fact that the discrepancy is strongest at the 
lowest Q 2 value, at higher Q 2 the resonance form factor reduces such contributions. 

Given this discrepancy in shape between the VGL calculations and the a\, data, the questions are: 1) how to 
determine the value of from the measured longitudinal cross sections <tl, and 2) what is the associated 'model 
uncertainty' in doing so? The difficulty is that there is no theoretical guidance for the assumed interfering background. 
This applies even if one assumes that the background is due to resonances: virtually nothing is known about the L/T 
character of resonances at W — 1.95 GeV, let alone their influence on <tl (via interference with the VGL amplitude). 

A. Summary of our Previous Extraction Method 

In our previous analysis |l3|, the following procedure was adopted. When fitting the value of A 2 (and hence F n ) 
for the separate i-bins the value of A 2 increases when —t decreases, since the data are steeper in \t\ than the VGL 
calculations. The value of extracted from the lowest |f |-bin, which is closest to the pole, was thus taken as a lower 
limit. 

An upper limit for F v was obtained by assuming that the background effectively yields a constant negative con- 
tribution to (7l- This background and the value of A 2 were then fitted together, assuming that the background is 
constant with t. The fitted contribution of the background was found to drop strongly as Q 2 increased from 0.6 to 
1.6 GeV 2 . Since in or this 'missing background' (i.e. the difference between the data and VGL model) decreases, at 
least for Q 2 =0.6 GeV 2 , with decreasing —t (Fig. [3|), and assuming that this also holds for ox, these assumptions give 
an upper estimate for F n . The best estimate for F^ was then taken as the average of the two results and one half of 
the average of the (relative) differences was taken as the 'model uncertainty'. 

However, the asumptions made in this procedure may be questioned. Firstly, the value of F v extracted from the 
lowest |i|-bin does not have to be a lower limit, and secondly, the assumption of a negative interfering ^-independent 
cross section in the upper limit calculation requires a special magnitude and phase for the interfering amplitude with 
respect to the VGL amplitude. 

B. Another Form Factor Extraction Method 

Since the publication of those results, we have looked at the discrepancy between the i-dependence predicted by 
the VGL model and the data in more detail by assuming, besides the VGL amplitude, a ^-independent interfering 
background amplitude, and fitting the latter together with the value of A 2 . The fitting uncertainty in A 2 varies 
between 5% and 18%, while the magnitude and phase of the fit background amplitude are very poorly constrained 
(uncertainties in the hundreds of percent). 

Although the fitting uncertainties are very large, the results of this exercise suggest an interfering amplitude whose 
magnitude decreases monotonically with increasing Q 2 , but whose phase with respect to the VGL amplitude does 
not necessarily result in a net negative cross section contribution to <tl, as has been assumed in the previous analysis. 
However, here also a special assumption was used, i.e. an interfering amplitude with a magnitude and phase that do 
not depend on t. Thus, determining F„ in this way is not a viable method, either. 

C. Adopted Form Factor Extraction Method 

Given that no information is available on the background, we are forced to make some assumptions in extracting 
F n from these data. Our guiding principle is to minimize these assumptions to the greatest extent possible. The 
form factor extraction method that we have adopted relies on the single assumption, that the contribution of the 
background is smallest at the kinematic endpoint t m i n . 

Our best estimate for F„ is thus determined in the following manner. Using the value of A 2 as a free parameter, 
the VGL model was fitted to each i-bin separately, yielding A 2 (Q 2 , W, t) values as shown in Fig. O The values of A 2 
tend to decrease as —t increases, presumably because of an interfering background not included in the model. Since 
the pole cross section containing F v increases strongly with decreasing —t, and the background presumably remains 
approximately constant, as suggested by the difference between data and VGL calculations for ctt, we assume that 
the effect of this background will be smallest at the smallest value of \t\ allowed by the experimental kinematics, 
|£ m in|- Thus, an extrapolation of A 2 to this physical limit is used to obtain our best estimate of F^. The value of A 2 
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FIG. 5: (Color online) Values of A 2 determined from the fit of the VGL model to each i-bin, and linear fit to same. The error 
bars reflect the statistical and f-uncorrelated systematic uncertainties. The additional overall systematic uncertainties, which 
are applied after the fit, are not shown. 



Q 2 W 
(GeV 2 ) (GeV) 


A 2 
(GeV 2 ) 




0.60 1.95 
0.75 1.95 
1.00 1.95 
1.60 1.95 


0.458 ± O.OSlto; 2 ,^ 
0.388 ± 0.038±g;ol| 
0.454 ± 0.034±g;g^g 
0.485 ± 0.0381^27 


0.433 ± 0.0171^36 
0.341 ± 0.022lg;^? 
0.312 ±0.016lo;o?g 
0.233 ± 0.014±g;g^ 


0.70 2.19 


0.627 ± 0.0581°;°^ 


0.473 ± 0.023±g;g|f 



TABLE II: A 2 and F n values from this work, and the reanalyzed data from Ref. [12[ using the same method. The first error 
includes all experimental and analysis uncertainties, and the second error is the 'model uncertainty' as described in the text. 



at t m i n is obtained by a linear fit to the data in Fig. [5) The resulting Aj and values are listed in Table HT1 The 
first uncertainty given represents both the experimental and the linear fit extrapolation uncertainties. 

The F-K values listed in Table [IT] correspond to the true values within the context of the VGL model if, and only if, 
the background vanishes at \t — t m i n \ = 0. Because of the uncertainty inherent in this assumption, we also estimate 
a 'model uncertainty' to account for the possible influence of the missing ingredient in the VGL model (background) 
at \t — t m i n \ = 0. Lacking a model for the background, we can only try to make a fair estimate of this uncertainty. 
This was done by looking at the variation in the fitted values of Aj when using two different assumptions for the 
background. We used the two cases considered earlier in this paper when trying to determine F n . The first case 
assumes the ^-independent negative background in addition to the VGL model used in Sec. IIII Al The second case 
assumes the interfering background amplitude with a ^-independent magnitude and phase discussed in Sec. IIII Bl 
However, here they are not used to determine F n , but only to estimate the model uncertainty in our best value of F^ 
determined above. 

The estimated model uncertainty is determined from the spread of the A 2 values at each Q 2 given by these two 
methods. Each effectively represents a different background interference with the VGL model. To keep the number 
of degrees of freedom the same in all cases, the background was fixed to the minimum \ 2 value determined in each 
of the above two studies, and A^ and its uncertainty was then determined in a one-parameter fit of the VGL model 
plus background to the ctl data. Since there is a strong statistical overlap between the two fits, the statistical plus 
random uncertainties of the data were quadratically removed from the A \ uncertainties. The model uncertainty at 
each Q 2 is then assigned to be the range plus fitting uncertainty given by the two methods, relative to the value of 
A 2 at t m i n . The resulting (asymmetric) uncertainties are listed as the second uncertainty in Table ITT1 
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The model uncertainty in the F„ value drops from about 20% at Q 2 =0.6 GeV 2 to about 5% at 1.6 GeV 2 . This is 
consistent with the fact that the discrepancy with the i-dependence of the VGL calculation is smaller for the larger 
values of Q 2 . It is also at least compatible with the idea that resonance contributions, which presumably have a 
form factor that drops fast with Q 2 , are responsible. The corresponding model uncertainty quoted in Ref. [l3| was 
approximately 8.5% at all Q 2 , but that was based on a more restrictive assumption on the background (in essence 
case 1). 

IV. RESULTS AND DISCUSSION 

Because of the arguments given above, the values presented in Table |TT] and Fig. [5] are our final estimate of 
from these data using this model. However, we stress that the primordial results of our experiment are the ctl cross 
sections. When improved models for the 1 H(e, e'-K + )n reaction become available, other (better) values of F^ may be 
extracted from the same cross sections. 

The present values for F n are between 7% and 16% smaller than our previously published values jl3| . which is 
about the combined experimental and model uncertainty. The largest difference is at Q 2 = 0.75 GeV 2 . On average, 
one quarter of the difference is because the final values of ctl are smaller than those of Ref. [l3| (see Fig. 0] for 
a representative comparison), and the remainin g th ree quarters are due to the F„ extraction method, the present 
method being closer to the method used in Ref. 13] to obtain the lower limit. 

0.5 H 1 1 1 1 1 h 




Q 2 (GeV 2 ) 

FIG. 6: (Color online) Q 2 F n data from this work, compared to previously published data. The Brauel et al. [T2| point has been 
reanalyzed using the F n extraction method of this work. The outer error bars for this work and the reanalyzed Brauel et al. 
data include all experimental and model uncertainties, added in quadrature, while the inner error bars reflect the experimental 
uncertainties only. Also shown are the light front quark model [y] (dash-dot), Dyson-Schwinger 0] (solid), QCD sum-rule [24[ 
(dot), dispersion relation [25[ (long-dash), and quark-hadron duality [2(J (short-dash), calculations. 

Analyses of other data at higher W indicate that the discrepancy with the ^-dependence of the VGL calculation is 
much smaller at higher values of W. The data from Brauel et al. [12|], taken at <3 2 =0.70 GeV 2 and a value of VF=2.19 
GeV, were reanalyzed using the F v extraction method presented here. The result is 0.4% higher than that obtained 
using the F n extraction method of Ref. 13]. This indicates that our F v extraction methods are robust when the 
background contribution is small, as appears to be the case at this higher value of W. 

The data from our second experiment [22j at IF~=2.22 GeV and Q 2 =1.6, 2.45 GeV 2 are also shown in Fig. [6] 
There, the VGL model adequately describes the ^-dependence of the <tl data, again indicating that the background 
contributions for <tl are smaller at higher W, even though the model strongly underpredicts the magnitude of ctt- In 
that case, the VGL model was fit to the full i-range of the ox data with only small fitting uncertainties. It is seen 
in the figure that the revised Q 2 =1.6 GeV 2 result from this work agrees well with that from our second experiment, 
taken at higher W and 30% closer to the ir + pole. The excellent agreement between these two results, despite their 
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significantly different t m i n values, indicates that the uncertainties due to the ir + electroproduction reaction mechanism 
seem to be under control, at least in this (J 2 -range. 

Fig. [6] compares our final data from this work and from our second experiment [22} to several QCD-based calcu- 
lations. The combined data sets are consistent with a variety of models. Up to Q 2 =1.5 GeV 2 , the Dyson-Schwinger 
calculation of Ref. [3], the light front quark model calculation of Ref. [6], and the QCD sum-rule calculation of 
Ref. are nearly identical, and are all very close to the monopole form factor constrained by the measured pion 
charge radius [?|. Such a form factor reflects non-perturbative physics. Our revised data are below the monopole 
curve. A significant deviation would indicate the increased role of perturbative components at moderate Q 2 , which 
provide in that region a value of Q 2 F^ «0. 15-0.20 only (2j. The dispersion relation calculation of Geshkenbein et al. 
[25| is closer to our results in the Q 2 =0.6-1.6 GeV 2 region, while still describing the low Q 2 data used for determining 
the pion charge radius. The quark-hadron duality calculation by Melnitchouk [26[ is not expected to be valid below 
Q 2 =2.0 GeV 2 . This is reflected in its significant deviation from the monopole curve at low Q 2 . To better distinguish 
between these different models, it is clear that especially higher Q 2 data, as well as more data at higher values of W 
in the Q 2 =0.6-1.6 GeV 2 region, are needed. Plans are underway to address both of these at JLab. 

V. SUMMARY AND CONCLUSIONS 

To summarize, the data analysis for our 1 H(e, e'ii + )n experiment at Q 2 =0.6-1.6 GeV 2 , centered at W = 1.95 GeV, 
has been repeated with careful inspection of all steps. The final unseparated cross sections presented here are in most 
cases consistent with our previous analysis within experimental uncertainties. After the magnifying effect of the L/T 
separation, the resulting ctt values arc slightly larger than before, and the ox values are correspondingly smaller. 
The experimental systematic uncertainties were critically reviewed, and are slightly smaller compared to the previous 
analysis. 



As before, we use a fit of the Regge model of Ref. [23[ to our <tl data to extract F„. The data display a steeper t- 
dependcncc than the model, which we attribute to the presence of longitudinal background contributions not included 
in the model. After revisiting our prior assumptions used to extract F„ from ol with the model, we conclude that 
some of our prior analysis assumptions were unwarranted. Therefore, we employ a revised F v extraction method 
which relies only on the assumption that the background contributions are minimal at t m j n . The resulting values are 
our best estimate of F^ from these data with this model, and are between 8% and 16% smaller than before, primarily 
due to the different extraction method. The Brauel et al. [l2j data at similar Q 2 but higher W are robust against 
our fitting assumptions, consistent with our expectation that a longitudinal background contribution not included in 
the Regge model is the cause of the discrepancy. 

The new analysis, in addition to providing our final F„ results for the Q 2 =0.6-1.6 GeV 2 range, gives an indication of 
the contribution of the analysis assumptions to the F v determination. A detailed analysis yields model uncertainties 
that decrease with increasing Q 2 and W. They are consistent with the differences in the values of F„ determined 
using the previous and present extraction methods. They indicate that, given the present electroproduction model, 
the uncertainty in the dtermination of F„ in this W and Q 2 range is of the order of 10%. 

The revised data indicate that for Q 2 > 0.5 GeV 2 , F^ starts to fall below the monopole curve that describes the 
low Q 2 elastic scattering data. These results are consistent with those of our second, more precise experiment at 
higher Q 2 and W [Hj]. Tne tw0 sets 01 data at Q 2 = L6 Gey2 

are taken with significantly different t m i„ , and so if the 
various form factor extraction issues were not being handled well by the VGL model, a significant discrepancy would 
have been expected to result. Their good agreement lends further credibility to the analysis presented here. It will 
be useful to acquire additional electroproduction data in the 0.5 < Q 2 < 1.5 GeV 2 range at higher W in order to be 
able to extract more precise values of F^ without the difficulties encountered here. 
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